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Absolute calibration of a wideband antenna and 
spectrometer for sky noise spectral index measurements 
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Alan E. E. Rogers, and Judd D. Bowman, 

A new method of absolute calibration of sky noise temperature using a 3-position 
switched spectrometer, measurements of antenna and low noise amplifier impedance with 
a vector network analyzer, and ancillary measurements of the amplifier noise waves is 
described. The details of the method and its application to accurate wideband 
measurements of the spectral index of the sky noise are described and compared with 
other methods. 



1. Introduction 

Accurately calibrated sky noise observations are 
needed to measure the spectral index of the diffuse 
radio background which is important in determin- 
ing the emission mechanisms in the Galactic and ex- 
tragalactic background. Theoretical models of the 
early universe predict small deviations in the spec- 
trum of the radio sky due to emission and absorption 
of neutral hydrogen at redshifts which place the 21- 
cm line in the 50 - 200 MHz range. In addition to 
the applications in radio astronomy accurately cali- 
brated sky noise measurements are used to measure 
the absorption of radio signals in the ionosphere us- 
ing riometers [Hargreaves and Detrick, 2002]. Also 
accurate calibration has applications in antenna pat- 
tern measurements and spectrum monitoring. The 
measurement of the antenna temperature is affected 
by the antenna impedance as any mismatch to the 
receiver results in a fraction of the sky noise power, 
T s k y , being reflected back into the sky so that the 
measured noise temperature from the antenna, T ant , 
is 



T an t — T s ky{l — |T| 2 ) 



(1) 



where T is the reflection coefficient. In addition out- 
going noise from the receiver is reflected back to the 
receiver by the mismatch. The "mismatch" can be 



corrected by using a mechanical tuner to adjust the 
antenna impedance to obtain the same impedance 
as the calibrated noise source used to calibrate the 
receiver. Under this condition any power lost due 
to the Low Noise Amplifier (LNA) having a differ- 
ent impedance is the same for the antenna as it is 
for the noise source. Several accurately calibrated 
measurements of the sky noise spectral index using 
this method were published in the 1960s and 1970s 
[Pauliny-Toth and Shakeshaft, 1962; Turtle et al, 
1962]. However this method is limited to a nar- 
row frequency range and requires constructing scaled 
antennas to cover a wide frequency range. With 
the development of broadband antennas, amplifiers 
and spectrometers new methods are needed for au- 
tomated measurements over wide frequency range. 

2. 50 ohm measurements at a fixed 
reference plane 

A vector network analyzer (VNA) provides an ac- 
curate measurement of the reflection coefficient ref- 
erenced to the standard 50 ohm impedance and con- 
sequently it is convenient to use a VNA to measure 
both the antenna and receiver input reflection co- 
efficient at the fixed reference plane defined by the 
50 ohm connection between the receiver and the an- 
tenna. In this case, with some algebraic manipula- 
tion equation (1) becomes 
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T sky (l-\T\')=T sky (l-\T a n\F\ 2 (2 



where 



{Z a + Zi) 



(3) 
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F = 



2)1/2 



(i-N 2 ) 
(i - r a r,) 



r„ = 



(Zg ~ 50) 

(Z a + 50) 
(Zi - 50) 



(4) 



(5) 



(6) 



[Z x + 50) 

Z a ,Zi and r a ,T/ are the impedance and reflection 
coefficients of the antenna and receiver respectively. 
Part of the complex factor F can also be viewed as 
the sum of the noise waves back and forth from the 
antenna to the receiver since the sum is a polyloga- 
rithmic series which converges when |r o rj| < 1 and 
has an exact solution 



where 



g(G(T amb + T C ai) + T ) 



G=l-N 2 



(11) 



(12) 



and g is the receiver gain, Tamb is the ambient tem- 
perature, T ca i is the excess noise temperature of the 
calibration noise. The calibrated receiver output, 
T 3p , is 



'3p 



T C al{Pant Pload) 



amb 



{Peal Pload) 

= T sky (i- irai 2 )^ 2 g- 1 

+ T u \T a \ 2 \F fc- 1 

+ (T c cos(cf>) +T s sin((f>))\r a \\F\G- 1 



(13) 



E( r ° r ') 4 = 77 



i=0 



(i-r„r,) 



(7) 



We can add the receiver LNA noise waves to equa- 
tion (2) to obtain the noise temperature, T reci from 
the LNA 



|ra| 2 )|F| 2 



T rec — T s ky{^- 

+{T c cos(cf>) +T s sin((j>))\r a \\F 



T u |r a | 2 |F| 2 



To 



(8) 



where T u is the uncorrelated portion of the LNA 
noise reflected from the antenna, T c and T s are the 
cosine and sine components of the correlated noise 
from the LNA reflected by the antenna and To is the 
portion of the noise which is independent of the LNA 
input. (j> is the phase of T a F. Equation (8) is equiv- 
alent to the noise wave formulation of Meys [1978] 
if the "reference" impedance is that of the LNA in 
which case F = 1. 

3. Receiver calibration 

A three position switch is used to provide the re- 
ceiver calibration. This switch sequentially connects 
the LNA to the antenna, a 50 ohm ambient load and 
a 50 ohm calibrated noise source. The power in each 

position, Pant, Pload, Peal IS 



Pant — S^T' 



Pload = g(GT a mb + Tq) 



(9) 



(10) 



4. Measurement of LNA noise wave 
parameters 

The LNA noise wave parameters, T Ul T c and T s 
can be measured by connecting an open (or shorted) 
low loss cable to the calibrated 3-position switched 
spectrometer in place of the antenna. The cable acts 
just like an antenna looking at an isotropic sky with 
temperature equal to the cable's physical tempera- 
ture. The cable reflection coefficient 
approximately 



r c , should be 



T c = L c exp{—iu)T) (14) 

where t is the 2-way delay, to is the frequency in 
radians/s and L c is the cable loss factor 



L c = 10 



-Jc/10 



(15) 



where lc is the cable loss in dB. For highest accuracy 
the open cable's reflection coefficient should be mea- 
sured with the VNA since most cables are not exactly 
50 ohms and the reflection coefficient is more com- 
plex than the expression above. With the open cable 
connected the calibrated receiver output is given by 
equation (13) by substitution of T s k y and r a with 
T am b and r c . Since the LNA noise parameters change 
gradually with frequency the rapid change of phase, 
(f>, afforded by a cable of sufficient length allows the 
noise wave parameters to be separately estimated 
via a least squares fit to the functions r a | 2 |F| 2 G _1 , 
co S (0)|r aWFlG- 1 and ain(^)|r a ||F \G~\ 
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5. Corrections for antenna losses 

If the antenna is lossless the calibrated sky noise 
temperature can be obtained from equation (13) 

T sky = [Tsp-TvlTaflFfG- 1 

- (T c cos(0) +T s sin(cl ) ))\r a \\F\G- 1 } 

x [(l-I^HlFpG- 1 ]- 1 (16) 

noting that when T a = 0, T s k y = T 3p . 

If the antenna loss factor is L the corrected sky 
temperature, T cs k y , is 



6. Model for choke balun loss 

The "choke" balun shown in Figure 1 is modelled 
by the circuit shown in Figure 2. In this model twice 
the ferrite impedance appears in parallel across the 
balanced antenna terminals so the reflection coeffi- 
cient measured at the unbalanced port, T, is 

T = Lexp(iujT)T b (22) 

where Tb is the reflection coefficient of the antenna 
impedance Z a , in parallel with twice the choke 
impedance, Zf/2, 



T C sky ~ (T U sky — T am b(L — 1)) / L (17) 

where T us k y is the uncorrected sky temperature and 



L = ltr'/ 10 



(18) 



where I is the loss in dB. For a typical antenna the 
loss is made of several components which include 
ground loss, resistive loss, transmission line loss and 
balun loss. The ground loss factor is one minus the 
fraction of the antenna pattern which receives ther- 
mal radiation from the ground. The resistive loss is 
one minus the fraction of the real part of the antenna 
impedance that comes from the conductor resistance 
including the skin-effect. If a transmission line and 
balun are part of the antenna the loss can be incor- 
porated into a 50 ohm attenuation in equation (8) 
which becomes 



- Wi-|r a | 2 )|F| 2 L 

+ T amb ((l-L)L|r a | 2 + l-£) 

+ T u \T a \ 2 \F\ 2 L 2 

+ (T c cos(4>) + T s sin(cl>)\r a \\F\L + To (19) 



where 



F = 



(i-N 2 ) 1 / 2 



(i-r a iu) 

L is the loss factor of the cable or attenuator and 



(20) 



T a = T aa exp(ibJT) 



(21) 



where r is the 2-way cable delay and T aa is the an- 
tenna reflection coefficient without the cable or at- 
tenuator. 



Z b = {l/Zf + 1/ZaY 



(Z b - 50) 



(23) 



(24) 



{Z b + 50) 

Since the noise power from impedance Z is propor- 
tional to its temperature times the real part of the 
impedance the loss factor is 



L = 



Re(Z a )\Z f \ 



(Re(Z a )\Zf\ 2 + Re(Zf)\Z a \ 2 ) 



(25) 



7. Example of calibrated sky spectrum 
measurement from 55 to 110 MHz 

A short observation of the sky spectrum was made 
with a broadband dipole and a 3-position switched 
spectrometer at a remote site away from strong TV 
and RF radio signals to test the calibration method. 
The switched spectrometer made the measurements 
of the following: 

1] 50 ohm resistor heated to 100 K above ambient 
2] An open cable about 3m in length 
3] Sky noise at a remote site 
VNA measurements were made of the following: 
1] Antenna through the balun at the remote site 
2] Input of the LNA through the 3-position switch 
3] Ferrite balun impedance and back-to-back loss 
4] Open cable used to derive the LNA noise pa- 
rameters 

Measurements of the heated resistor were used to 
calibrate the internal noise diode in flatness and ab- 
solute scale. To remove any uncertainty in the tem- 
perature of this "hot" the measurement was checked 
using a HP346C noise source with calibration trace- 
able to NIST. 



4 
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The wideband dipole used for the test is based on 
the "Fourpoint" design of Suh et at [2004]. This 
antenna which was made of aluminum rod elements 
instead of panels is shown in Figure 3. Figure 4 
shows a block diagram of the system used to measure 
the sky noise spectrum which includes the antenna, 
balun, 3-position switch, LNA and spectrometer. It 
also shows a noise source used to inject added noise 
below 55 MHz which conditions the analog to digi- 
tal converter (ADC). This "out of band" condition- 
ing improves the linearity and dynamic range of the 
ADC. 

Figure 5 shows the measured reflection magnitude 
and phase along with a Fourier series least square fit 
to the data to reduce the noise in the VNA data and 
provide accurate interpolation between the measure- 
ment. 

Figure 6 shows the measured reflection coefficient 
of the LNA. In this case the value of fitting a Fourier 
series is clearly evident as the VNA measurements of 
the LNA were make with a spacing of 10 MHz. 

Figure 7 shows the spectrum of the calibrated 
3-position switched spectrometer connected to the 
open cable with 2-way delay of 29 ns and 0.7 dB 
loss. The solid curve is the spectrum derived from 
the best fit LNA noise parameters shown in Figure 
8. The frequency dependence of the noise parame- 
ters were constrained to a constant plus a slope. At 
80 MHz the values of T U ,T C ,T S were 35, 9, and 10 
respectively. 

Figures 9 and 10 show the calibrated sky spectrum 
for 2 sequential observations taken at the remote site 
in West Forks, Maine on 8 July 2011 at 17 UT. In 
the first the antenna was connected directly to the 
3-position switch and in the second the antenna was 
connected through a 6 dB attenuator. These cali- 
brated spectra were derived using equation (18) or 
equation (21), in the case of the use of the added at- 
tenuation, followed by equation (19) for the correc- 
tion of the antenna loss factor. The balun loss factor 
was calculated using equation (25) assuming a ferrite 
choke impedance of 400 ohms for which the circuit 
model gives a back-to-back loss of 0.52 dB for back- 
to-back baluns consistent with the measured value. 
A loss factor of 0.93 was assumed for the ground loss 
of 7% based on the results of an electromagnetic sim- 
ulation of the antenna and ground plane. The overall 
loss factor is the product of the balun and ground 
loss factors. The antenna was assumed to have no 
resistive loss. 



8. Determination of spectral index 

An estimate of the spectral index is obtained by 
finding the the function, 



T sky (f)^T skv (f0)(f/f0)- s P 



(26) 



where T s k y (f) is the sky temperature as a function 
of frequency /, sp is the spectral index and /0 is 
the "reference" frequency, that fits the calibrated 
sky noise spectrum using weighted least squares for 
which portions of the spectrum with added interfer- 
ence are given zero weight. The best fit spectral in- 
dices using the calibrated sky spectra shown in Fig- 
ures 9 and 10 were 2.50 +/- 0.05 and 2.53 +/- 0.05 re- 
spectively. These values are in good agreement with 
measurements made from 100 to 200 MHz by Rogers 
and Bowman [2008] in which the antenna was con- 
nected to the antenna via a long cable so that cor- 
related noise terms in equation (19) are fitted and 
removed. The disadvantage of this method is that a 
long cable has to be used to provide the rapid change 
of phase of r a which results in added loss. In addition 
the "ripples" in the measured spectrum are assumed 
to arise only the mismatch and cannot be separated 
from fine scale structure which might be present in 
the sky spectrum. In this method the antenna re- 
flection coefficient was measured by injecting added 
noise which was then reflected by the antenna and 
measured by the spectrometer. But the antenna re- 
flection coefficient could have been measured with a 
VNA. 

9. Sensitivity of calibrated sky temperature 
to VNA measurement errors 

The sensitivity to measurement errors in the ancil- 
lary measurement of antenna and LNA reflection co- 
efficient, balun ferrite impedance and ambient tem- 
perature was determined by using equation (19) to 
generate noiseless "reference" data assuming a sky 



Table 1. Sensitivity to changes in ancillary measurements 81 



Parameter changed 


rms (K) 


rms2 (K) 


r a amplitude 


1.8 


0.09 


r o phase 


0.25 


0.06 


T[ amplitude 


0.13 


0.06 


T; phase 


0.19 


0.04 


Zj, amplitude 


0.38 


0.09 


Zb phase 


0.004 


0.001 


T amb temperature 


1.2 


0.8 



rms2 is for antenna with -20 dB reflection coefficient 
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temperature of 1500 K at 80 MHz and a spectral 
index of -2.5. New data was then generated by per- 
turbing the ancillary VNA data, hrst by 1% in am- 
plitude, and then by 0.01 radians in phase. A test of 
the sensitivity to the ambient temperature was also 
made by changing the assumed ambient temperature 
by 1 K. These two data sets were then processed to 
obtain the calibrated sky temperature and after re- 
moving the best fit spectral index their rms differ- 
ences computed. The results with an attenuator of 
10 dB placed between the antenna and the 3-position 
switch are given in Table 1. 

For long integrations for which the systematics 
still dominate an attenuation of about 10 dB is opti- 
mum for the high sky noise temperature at 80 MHz. 
Using more attenuation makes little difference to the 
sensitivity to errors in the ancillary measurements 
and using no attenuator increases the sensitivity to 
error in the reflection coefficients by about a factor 
of 2. It is noted that the most critical measurement 
is that of the antenna impedance. However the sen- 
sitivity to antenna impedance can be reduced if an 
antenna with a better match over the wide frequency 
range can obtained. For example a large reduction in 
sensitivity could be obtained with an antenna whose 
the reflection coefficient is below -20 dB from 55 to 
110 MHz as shown by the numbers in the last col- 
umn of Table 1. A reduction of the sensitivity to 
antenna reflection coefficient by a factor of about 3 
can be obtained with the "Fourpoint" used in the test 
if the frequency range is limited to 60 to 90 MHz. We 
note that precision impedance analyzers are available 
which may achieve higher accuracy than a VNA in 
the frequency range below 110 MHz. 



10. Conclusion 

The accurate calibration method described has the 
potential to make wideband sky noise spectral mea- 
surements with greatly reduced level of systematics. 
Simulations of performance look encouraging but em- 
phasize the very high accuracy needed in the mea- 
surements of reflection coefficient. 
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balanced input 




Figure 1. Ferrite choke balun 
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Balun circuit model 



Figure 2. Circuit model of ferrite choke balun 
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Aluminum 




Figure 3. Fourpoint antenna 0.43A tip-to-tip dipole 0.2A above ground at 55 MHz 
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Figure 4. Block diagram of complete system for sky noise measurements 
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Figure 5. Antenna reflection coefficient 
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Figure 6. LNA reflection coefficient 
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Figure 7. Spectrum of open cable 
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Figure 8. LNA noise parameters. Top curve is the 
phase of the correlated noise. The middle curve is uncor- 
rected noise and the bottom curve is the magnitude of 
the correlated noise 
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Figure 9. Calibrated sky noise spectrum without attenuation 
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Figure 10. Calibrated sky noise spectrum with 6 dB attenuation 



